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INTRODUCTION 

Present U.S. Air Force [1] and proposed U.S. civil [2] airworthiness 
regulations are based on considerations of "damage tolerance" in aircraft 
structures. The underlying philosophy for these regulations is to 
acknowledge that accidental or normal service- induced damage is inevitable 
and that periodic inspections are required to detect such damage. Air- 
worthiness is then assured by demonstrating that damage that escapes one 
Inspection will not grow to critical size before the next inspection. Two 
evaluations must be made: first, the rate of damage propagation under 

expected service loading; and second, the residual static strength in the 
presence of significant damage. Both of these evaluations employ fracture 
mechanics analyses and predictions. 

Unfortunately, most fracture mechanics developments have concentrated 
on analysis and prediction of behavior in brittle materials and relatively 
simple heavy sections for which simple linear elastic analysis are adequate. 
In contrast, aircraft structures usually employ somewhat ductile materials, 
thin gages, and complex configurations. Further, complex load histories 
and hostile environments are generally features of most aircraft design 
problems. Consequently, special considerations are required to treat these 
problems adequately. 

NASA-Langley Research Center has conducted a variety of research tasks 
to extend the capabilities of fracture mechanics to deal with some of 
these complexities. The purpose of this paper is to describe the current 
stage of development of these capabilities. 


LIMITATIONS OF LINEAR ELASTIC FRACTURE MECHANICS 
Fracture 

Griffith [3] is generally cited as the earliest investigator of fracture 
mechanics. He was concerned with the behavior of glass, an inherently 
brittle material. Irwin [4] expanded on Griffith's principles to deal 
with fracture problems in welded steel ships that failed unexpectedly and 
catastrophically at cold temperatures. Under the conditions prevailing in 
the ship application, the material behaved in what was regarded as a brittle 
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fashion. Both Investigators recognized that, although the analytical 
framework they applied was physically appropriate for linear-elastic, or 
brittle, behavior, the observed fracture strength was different by orders 
of magnitude from the values predicted by physical conditions alone. 

Thus, practical applications of fracture mechanics generally •'equire 
empirical correction factors to account for the actual behavior of the 
materials in question. During the last three decades, when fracture 
mechanics has developed into a recognized discipline, a frequently employed 
device Is to force brittle behavior in a given material by making research 
specimens thick enough to develop essentially a plane-strain stress state 
at the crack tip. This device makes the linear elastic theory more nearly 
capable of correlating results. Unfortunately, many practical materials 
require that test specimens be up to 200 nr* thick to produce "valid" test 
results. At least for aircraft structures, the practical range of thick- 
nesses Is about 1 to 35 mm. Thus, linear elastic fracture mechanics cannot 
deal directly with many practical problems. 

To correct this deficiency, a wide variety of fracture analysis tech- 
niques have been developed. Several of these are described in other papers 
In the present symposium. In most cases a quasi-rigorous rationale is 
developed to account for some recognized phenomenon that makes real 
behavior different from the linear elastic behavior, and usually an 
empirical adjustment is ultimately required to fit the equations to actual 
data. The present paper follows this same pattern. 

One distinguishing feature of the NASA work is best explained with the 
aid of Fig. 1. Ideal brittle behavior is depicted in Fig. 1(a). A center- 
crack tension specimen contains a crack of length Cq and is subjected to 
monotonically increasing load P (plotted as the ordinate). Because all 
behavior is linear, the load may be increased without appreciable crack 
growth until a critical value Perit reached, at which time the speciiren 
suddenly fails. As indicated earlier, this behavior is rare, even in 
materials regarded as brittle. 
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a. Brittle b. Practical 

Fig. 1 Fracture of ideal and practical materials 
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Practical materials behave as depicted In Fig. 1(b). A similar 
specimen with the same size original crack Is subjected to steadily 
increasing load. At point 1, the crack begins to grow slowly and stably. 
Depending on the material and configuration, considerable crack growth and 
local plasticity accon^any Increasing load until an instability is reached 
at point 2. The various "Improved" fracture analysis methods are directed 
toward answering an assortment of basically different questions. Some 
characterize the conditions for initiation of crack growth (point 1), 
others characterize the instability point (2), and do so in conjunction with 
identifying the corresponding crack length at the instant of instability, 
still others define the complete curve between 1 and 2 . 

The NASA work seeks to characterize the maximum load supported by the 
specimen just before failure and to associate that load with the initial 
crack length. This relationship is. of course, not located on the crack 
growth curve, but it answers a question of fundamental significance to 
designers and operators of practical hardware. 


Fatigue Cract Propagation 

Paris [5] was one of the first to apply fracture mechanics principles to 
the correlation of crack propagation data for tests under constant ampli- 
tude loading. He employed linear elastic analyses and proposed a simple 
power law to relate the crack growth per cycle with the range of the stress 
intensity caused by the load fluctuation in a given cycle for the crack 
configuration in question. Since that time, the simple power law has 
been shown to require adjustments to account for: an apparent threshold 
below which cracks appear not to grow; stress intensities approaching the 
critical intensity where unstable crack growth is expected; commonly 
observed effects of n»an and maximum stresses; and nonlinear effects 
observed when a complex time history of loads is employed. 

Investigations to account for these pherotrena employ a wide variety 
of empirical or semierripirical methods with varying degrees of success. 

As is explained in a later section, the NASA work is based on consideration 
of crack closure during the load cycle, a basic phenomenon first observed 
and quantified by Elber [6]. 


Effects of Complex Configuration 

All analyses of crack behavior require a basic stress analysis of the 
specimen or structural member under consideration. Generally, the local 
stress intensity is derived as a function of the shape of the part, the 
loading, the mode of crack deformation, dimensions of the crack (length and 
depth) relative to dimensions of the part (width and thickness), and 
plane strain or plane stress conditions at the crack front. Several 
catalogs of such solutions have been published [7-9] and are widely used. 
However, some published solutions disagree significantly or fail to 
account for some significant, but mathematically intractable, feature. 
Recent developments of high-speed computers have led to solutions for much 
more complicated configurations, including several in three dimensions. 
Although NASA work has contributed to the analysis of several unique cases 
[10-13], these will not be discussed in this paper. However, a later 
section deals with the analysis of complex structural assemblies to illus- 
trate how such configurations can be treated. 
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A TWO-PARAMETER FRACTURE CRITERION 

As Indicated earlier, NASA work on fracture has sought to predict the 
maximal load that can be carried by a cracked part. Because linear elastic 
analyses, by themselves, were Incapable of dealing with fracture in 
practical cases where considerable stable crack growth and plasticity 
occurred, a new criterion was developed [14-15]. Two parameters were 
Involved In the final relation, which was derived from a combination of 
stress Intensity and stress redistribution considerations. The criterion 
Is expressed as 





when 


< a 


ys 


( 1 ) 


where 


Kjg Is the elastic stress Intensity at failure. 

Is the nominal stress on the net section at failure, 
n 

S Is the nominal stress required to produce a plastic hinge 
** on the complete net section, using the ultimate tensile 
strength of the material as the maximum stress. 

Kp Is an empirically adjusted material constant having the 
sanK dimensions as an elastic stress Intensity. 

n Is an empirically adjusted material constant whose value 
nwst He between zero and one. 


Because two constants require adjustment, at least two (preferably 
considerably different) tests must be conducted to characterize the 
material. The parameter m Is multipled by the stress ratio Sn/S^ and is 
thought of as representing a plasticity correction. Its value is 0 if the 
material Is essentially brittle and 1 if the material is highly ductile. 

The parameter Kp Is thought to be an inherent fracture strength, dependent 
on thickness, but independent of in-plane configuration and loading. 

This criterion has been applied in the analysis of a large quantity of 
data taken from the published literature and from NASA tf*sts. A partial 
list of parameters represented In these data follows: 

Crack configurations; 

Center crack tension 
Contact 
Notch bend 
Surface cracks 
Corner cracks at holes 

Through-cracks and surface cracks In pressurized cylinders 

Materials: 

Steels Brasses 

Titanium alloys Magnesium alloys 

Aluminum alloys 

Thickness range: 0.5 mm to 50 mm 

Width range: 25 mm to 1250 mm 
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In all cases, the criterion depicted the appropriate trends among 
related sets of data. When two or more sets of data were available for the 
same material and the same thickness but syst&natically different specimen 
configurations, the criterion predicted failure loads for one set when the 
fracture parameters obtained for another set were used. Two representative 
treatments [16-17] are shown In Figs. Z and 3. 


t • 3.2 mm 



a. Curves fitted to data b. Curves predicted from 

other data 


Fig. 2 Correlation of fracture tests of Hiduminium H-48 alloy 



Fig. 3 



(b) Curve predicted from 
other data 


Correlation of fracture tests of 2219-T87 aluminum alloy 
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Figure 2 Is for two series of tests on a British aluminum alloy, 
called Hidufflinum H-48 [18]. Center cracked tension specimens with 
several widths and various crack lengths were tested. The two paran«ters 
were determined for best fit with all the results shown in Fig. 2(a). 

These same values of the fracture paranwters were then employed to predict 
the behavior of the compact specimen results shown in Fig. 2(b). In Fig. 3 
the results of tests of 2219-T87 aluminum alloy specimens with surface and 
comer cracks [19] are treated in similar fashion. In Fig. 3(a) the 
two parameters were determined to achieve a best fit with test results from 
surface-cracked tension specimens. Figure 3(b) shows that predictions from 
the criterion and the parameters found in Fig. 3(a) agreed well with 
results of tension tests on specin»ns having corner cracks near holes. 

Although no physical reason can be cited for establishing a correla- 
tion between the values of Kp and m, a cross-plot of these two parameters 
for a wide variety of materials and specimen configurations is shown 
In Fig. 4. The values of Kp were normalized by E, the elastic modulus for 
each material, to permit comparison of several classes of alloys. Although 
no generality of this relationship is claimed at this date, the points fall 

/AKp\ 

Into a reasonably small scatter band. The simple relation m = Tanhi-r-j 
has been plotted to represent this correlation. ^ ^ ' 



O ALUMINUM 
□ STEEL 
A TITANIUM 


SURFACE 

CRACKS 


• ALUMINUM 
■ STEEL 
▲ TITANIUM 


THROUGH 

CRACKS 


Fig. 4 Apparent relation between Kp and m 


FATIGUE CRACK PROPAGATION 
Constant-Amplitude Loading 

As indicated earlier, the direct application of the simple power law to 
predict fatigue crack propagation is inadequate for several reasons. Elber 
[6] was the first to identify and quantify the phenomenon of crack 
closure. He showed persuasively that, in materials that display appre- 
ciable ductility, the crack surfaces near the crack tip may close before a 
specimen is completely unloaded. Such behavior is not expected in elastic 
specimens until at least an Infinitesimal compression load is applied. 

The reason for the crack closure is depicted in Fig. 5. The presence 
of a plastic zone ahead of a crack tip (shown with double cross hatch) is 
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generally accepted as a feature of crack behavior in practical materials. 
The feature that had been comnnnly ignored was that the material adjacent 
to the crack surfaces near a crack tip had once been in this plastic zone, 
and thus, had been defomttd in tension to a size larger than the original 
size. Consequently, as the load on a specimen is reduced, this material 
forces the crack surfaces to meet, to transfer compressive forces across 
the Interface, and perhaps to produce reversed plasticity (compression) in 
the previously deformed material, all before the load is reduced to zero. 
Many investigators have conducted detailed experinffints tc convince them- 
selves and ethers that this phenomenon exists. The intuitive consequence 
of this behavior is that those portions of a load cycle during which the 
crack tip is closed will not contribute the crack propagation. 



In his early work, Elber [20] measured crack opening displacements 
near the crack tip to quantify the "crack opening stress level." He found 
a reesonably consistent agreement between test data on aluminum alloy 
specinen and the relation; 

So * Smax(o-5 + O.IR + O.AR^) if -0.1 < R < 0.7 (2) 

Where 

Sq is the stress level at which the crack tip opens 
Sfliax maximum stress applied 

R is the ratio S,nax 
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More recently, Newman [21] analyzed this behavior using a two-dimensional, 
elastic-plastic, f1n1te-el«nent analysis and showed quantitatively that 
cracks open and close at predictable stress levels. The finite-element 
analysis was shown to be quite accurate but was very con^Hcated and 
required a large computing facility. More recently, Newman [22] analyzed 
this behavior with a simple strip-yield model that requires only a small 
con^uter and can acconmwdate many hundreds of thousands of variable load 
cycles. The region near the crack tip was modeled by bar elements having 
dimensions as suggested by the small rectangles In Fig. 6. An elastlc- 
perfectly-plastlc stress-strain curve was assumed. The yield stress 
level was adjusted to a value 15 percent greater than the actual yield 
stress of the material to compensate for some biaxial Ity that Is expected 
In the highly stressed region. The rectangles In Fig. 6 show the. .extent 
of plastic deformations during the present and preceding load cycles. 

The vertical gaps between rectangles depict the crack opening. In the 
analysis, the amount of crack growth 1n a given cycle Is introduced Into 
the calculation by the analyst and is based on the effective stress 
range (Eq. (2)) and the simple power law for crack growth In the material. 
The stress distribution In the lower left of the figure shows the local 
stress to be zero where the crack Is open, equal to the effective yield 
stress In the plastic zone, and the elastic values elsewhere. 

Upon unloading (right side of Fig. 6), the release of elastic strains 
throughout the model forces the crack surfaces to approach each other, 
touch, transfer compression force and undergo reversed plastic deformation 
when appropriate. The stress distribution in the lower right shows the 
local stress to be zero where the crack Is open, to equal the effective 
compressive yield stress where reversed plasticity has occurred (u) and to 
have elastic values elsewhere. This nwdel also accounts for compressive 
yielding at the crack surfaces, If appropriate. 



Fig. 6 Calculations of local stresses and deformations 
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This model has been used to predict the craci' -opening stress level for 
any con4>1nat1on of maximum and mlnlimim stress and /or prior cyclic load 
history. Figure 7 shws the crack-opening stresses nonnallzed by 
plotted against the ratio of to the effective flow stress (cq). The 
effective flow stress Is asstmied to be the average between yield stress and 
ultimate tensile strength. 

EXPERIMENTAL lEQ. 2) 

• FINITE-ELEMENT ANALYSIS [REF. 21j 

STRIP -YIELD MODEL [REF. 22] 




Fig. 7 Experimental and predicted crack opening stresses 


The results In Fig. 7 were used to analyze fatigue crack propagation 
data for 7075- T6 alumiman alloy center-cracked sheet specimens [23] 
with a broad range of R values and maximum stresses. In Fig. 8 the basic 
data are shown plotted in the conventional form of growth per cycle against 
the range of the elastic stress intensity factor. The several discrete 
lines of test data reflect the R effect. 

The same data are shown in Fig. 9 except that the range of stress 
Intensity was adjusted for two effects. First, crack closure effects were 
accounted for by; 

A*^ff “ (^x ‘ So)\/ifS F (3) 

Where 

Sq Is the crack opening stress level 

F Is the usual paran«ter that accounts for specimen configuration. 
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1 10 100 
AK. MPa 


Fig. 8 Crack growth data from tests at various R values, 
7075>T6 aluminum alloy sheet 



1 10 ion 

^ MPa 


Fig. 9 Correlation of crack growth data from tests at 
various R values, 7075-T6 aluminum alloy sheet 
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Second, the effect of high stresses was accounted for by*. 



Is the stress Intensity corresponding to Sma^ 

Kp Is the fracture parameter discussed earlier. 

This equation has a form similar to an empirical equation first proposed by 
Forman [24] to account for mean stress and high stress effects. 

All the data are seen to fall in a very much narrower scatter band, 
Indicating that the adjustments ^or crack closure, R effects, and maximum* 
stress effects are appropriate. A single power law has been fitted to 
these transformed data. Obviously, this law must not be used below a 
threshold value, which was not established for these data. 

A similar treatment of data for 2024-T3 aluminum alloy sheet sped* 
mens [23] Is shown in Fig. 10 and illustrates a similarly excellent 
correlation. Few other sets of data exist that cover a wide enough range 
of stress ratios and stress levels to make a convincing argument for this 
analysis. However, the n«thod has been applied to several other smaller 
sets of data with encouraging results. 



Fig. 10 Correlation of crack growth data from tests at 
various R values, 2024-T3 aluminum alloy sheet 
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Var1abl6>Ainp11tude Loading 


Essentially all practical crack propagation evaluations require considera- 
tion of load histories with variable amplitude loading. A strong nonlinear 
behavior usually called a “delay" generally follows high loads in the time 
history. Figure 11 Illustrates a particular case. Center-cracked 2024-T3 
aluminum alloy sheet specimens [22j were subjected to the simple two- 
amplitude loading schedule shown in the Inset. The crack length is plotted 
against the nufrt>er of load cycles. The dashed curve at the left is the 
behavior observed ^f no spike load had been applied. The symbols labeled, 
"Test 1“ and "Test 2" represent results of two identical ttsts in which the 
spike load was applied at N 0. These two test results are not identical 
because of normal scatter in behavior. However, both these tests show a 
"delay" of 100, JOO or more cycles during which essentially no crack propa- 
gation was measured. 



rig. 11 Predicted and experimental crack growth in material 
subjected to variable amplitude loading, 2024-T3 aluminum 
aikOy sheet 


The curve labeled "Predicted" was calculated by the previously described 
procedure, except that the crack opening stress level associated with the 
spike load (computed from the simple strip-yield model) was used to define 
the effective range of the stress intensity. Because the crack opening 
stress was almost equal to the maximum stress in cycles following the spike 
load, essentially no crack propagation was predicted for the first 100,000 
cycles of the lower tress amplitude. The prediction agrees satisfactorily 
with the test data. 

A more complex, but practical, case of crack propagation under variable 
amplitude loading is shown in Fig. 12. For this example, a 2024-T3 alumi- 
n«mi alloy sheet specimen containing a center crack was subjected to the 
simulated service load history shown in the inset. For this example, the 
highest stress In the time history was 140 MPa. The stress on the specinten 
was reduced to zero periodically during each test to simulate the unloading 
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of a point on the lower surface of an aircraft wing once for each flight. 
Tte nut^r and magnitude of stresses In a given flight was varied to agree 
with statistical est1n»tes of loading frequencies experienced In service of 
a civil transport aircraft. 


ONE 



Fig. 12 Crack propagation under simulated flight loading, 
2024>T3 aluminum alloy 


The dashed curve on the left shows the behavior expected if crack 
propagation were to proceed linearly for each stress excursion (no load 
interaction). The test data show that the speciron survived almost 
70,000 simulated flights instead of the few thousand flights expected if no 
load interaction were present. The predicted crack growth was again based 
on crack closure predicted by the strip-yield model. Although the pre- 
dicted number of simulated flights was too high for crack growth between 1 
and 2 rmi, the predicted and experimental curves agree well for longer 
cracks. Further refinement of the analysis is expected to improve the 
agreement for short and intermediate crack lengths. 

Because a given gust or maneuver load on an aircraft wing may cause a 
variety of local stress levels depending on detailed configurations, this 
analysis v«s repeated with the highest stress in the tii.ie history syste- 
matically varied. The predicted and experimental numbers of flights to 
failure are plotted as bar graphs in Fig. 13. Three values of maximum 
stress v^re assumed: 140, 170, and 200 MPa. The agreement between predic- 
tions and experimental results was excellent in each case. 
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MAXIMUM STRESS LEVEL IN SPECTRUM 


Fig. 13 Crack propagation life under simulated 
flight loading, 2024-T3 aluminum alloy 


CRACK PROPAGATION AND FRACTURE IN BUILT-UP STRUCTURES 

Aircraft wing structures are inherently composed of complex assemblies of 
sheets, stiffeners, splices, and doublers in order to achieve maximum 
strength with minimum mass and enhance damage tolerance. Thus, effective 
design requires that the behavior of such a structure be analyzed to 
determine rates of crack propagation and residual strength under service 
loadings. 

Poe [25] has developed an analysis for calculating the stress 
intensity, rivet force, and stringer loads in a stiffened panel containing 
an initial crack. The sketch at the top of Fig. 14 depicts a typical con- 
figuration that has been analyzed. The analysis accounts appropriately for 
the several configuration parameters listed in the figure. The ordinate of 
the plot has two scales: (1) the stress intensity in the sheet normalized 

by the stress intensity for an unstifferied sheet containing a crack of 
equal length; and (2) the load in the central stringer normalized by the 
load in that stringer if no crack were present. The abscissa is the half- 
length of the crack and corresponds to the crack length shown in the sketch 
According to the analysis, the stress intensity in the sheet is reduced sig 
nificantly as the crack tip passes each stiffener. This illustrates that 
stringers can influence crack growth rates significantly. The amount by 
which the stress intensity is reduced depends on the specific values of 
each of the configurational parameters listed. 
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Fig. 14 Stress analysis of stiffened panels containing damage 


Predictions from this analysis have been compared with rates of crack 
propagation observed [26] in fatigue tests conducted on systematically 
varied configurations and stress levels. The trends in the data were 
always predicted appropriately and good quantitative agreen^nt was achieved 
for each case. 

The upper curve in Fig. 14 illustrates that the central stringer 
carries a progressively higher share of the load in the structure as the 

crack grows. This load levels off when the crack passes other stringers, 

because those stringers then carry higher loads. 

The residual static strength of a stiffened structure containing a 
crack has been assessed on the basis of this analysis. Such an assessment 
is described in Fig. 15. In this case, the stress required to cause static 

failure in the sheet or in the stringer is plotted as the ordinate. The 

two-parameter fracture criterion, described earlier, was used to predict 
the falling stress level in the 2024-T3 aluminum alloy sheet, and the ulti- 
mate tensile strength of the stringer material was used to predict the 
failing stress level for the stringers. Stringer 1 was assumed failed at 
the outset. The crack length in the sheet is plotted on the abscissa and 
alined with the sketch of the panel in the inset. 

If the stress is progressively increased on a panel containing the 
crack length shown, the predicted sequence of failure events is shown by 
the path marked with arrows. No failure occurs until the applied stress 
reaches the level of point A, when the sheet is expected to fail. Because 
the residual strength is less than A for all crack lengths between A and B, 
the failure will extend to B, but will be arrested there, because more 
stress must be applied for crack lengths longer than B. As the stress is 
increased the crack will grow until point C is reached. At C, stringer 2 
will have reached its ultimate tensile strength and will fail. Because the 
forces formerly carried by this stringer are suddenly transferred to the 
sheet, the load carrying capability of the sheet is suddenly reduced to the 
level labeled 0. Thus, unless the load is removed instantaneously, the 
crack will quickly extend to point E where stringer S3 will fail and all 
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othtr parts of the cross section will fail in turn, because the residual 
strength never again exceeds the level of C. 


S-1 S-2 S-3 



FAIL-SAFE 

STRENGTH 


The relative stiffnesses of the sheet and stringers, stringer spacings, 
rivet spacings, and other configurational parameters modify the relative 
positions of the several curves in Fig. 15. Thus, effective design requires 
that a range of these parameters be evaluated to produce the most damage 
resistant configuration with minimum weight and cost. 

To facilitate this analysis, a first-generation computer program [27] 
has been written to deal with this design problem. The program predicts 
both the rate of crack propagation and the residual strength of a built-up 
configuration under service loadings. From these considerations the best 
configuration may be chosen and the allowable stress levels may be 
adjusted to provide the desired characteristics. Other computer programs 
are available to deal with other disciplinary concerns, such as original 
static strength, stiffness, flutter, and buckling, so that the design satis- 
fying several criteria may be achieved with minimum weight penalty. 


CONCLUDING REMARKS 

This paper has reviewed systematic research that has led to improved 

1. Stress intensity solutions for laboratory specimens 

2. Fracture criterion for practical materials 

3. Crack propagation predictions that account for mean stress and high 
maximum stress effects 

4. Crack propagation predictions for variable-amplitude loading 

5. Prediction of crack growth and residual strength in built-up 
structure assemblies 
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These several capabilities have been Incorporated Into a first-generation 
computerized analysis that allows design for damage tolerance and trade- 
offs with other disciplines to produce efficient designs that meet current 
airmrthlness requirements. 
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